Squeezed-light generation in a nonlinear planar waveguide with a periodic corrugation 
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Two-mode nonlinear interaction (second-harmonic and second-subharmonic generation) in a pla- 
nar waveguide with a small periodic corrugation at the surface is studied. Scattering of the interact- 
ing fields on the corrugation leads to constructive interference that enhances the nonlinear process 
provided that all the interactions are phase matched. Conditions for the overall phase matching arc 
found. Compared with a perfectly quasi-phase-matched waveguide, better values of squeezing as 
well as higher intensities are reached under these conditions. Procedure for finding optimum values 
of parameters for squeezed-light generation is described. 

PACS numbers: 42.50.-p Quantum optics, 42.65.Ky Frequency conversion, 42.65.Wi Nonlinear waveguides, 
42.70.Qs Photonic bandgap materials 
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I. INTRODUCTION 

Since the pioneering work by Armstrong [lj on the pro- 
cess of second-harmonic generation has occurred, spatio- 
temporal properties of the nonlinearly interacting clas- 
sical fields have been studied in detail by many authors 
both theoretically and experimentally. A new impulse in 
these studies has occurred when people understood that 
this process can give rise to the fields with nonclassical 
properties (for a review, see, e.g. Q)- Namely light with 
electric-field amplitude fluctuations suppressed below the 
limit given by quantum mechanics can be generated both 
in the pump and second-subharmonic fields. Also light 
with nonclassical photon-number statistics can be ob- 
tained - pairwise character of photon-number statistics 
[H, 0| generated in the spontaneous process of second- 
subharmonic generation has been observed 

It has been shown that the best conditions for 
squeezed-light generation in homogeneous nonlinear me- 
dia occur, provided that the nonlinear two-mode interac- 
tion is perfectly phase matched. Under these conditions, 
the principal squeeze variance of the second-subharmonic 
field can asymptotically reach zero when the gain of the 
nonlinear interaction increases. On the other hand, the 
pump-field principal squeeze variance cannot be less than 
0.5 [fj. If large values of the nonlinear phase mismatch 
are allowed, this limit can be overcome due to a nonlinear 
phase modulation, as suggested in 0, However, the 
generated signal is very weak. 

The most common method how to compensate for 
the natural nonlinear phase mismatch that occurs in 
commonly used nonlinear materials is to introduce an 
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additional periodic modulation of the x susceptibil- 
ity using periodical poling [13, [llj. Several meth- 
ods for the additional modulation of the local ampli- 
tude of this quasi-phase-matched interaction have been 
developed [12| , These methods allow to reach a spec- 
trally broad-band two-mode interaction and so femtosec- 
ond pumping of the nonlinear process is possible [l3j . 

In order to effectively increase low values of the nonlin- 
ear interaction in real materials, configurations in which 
a nonlinear medium is put inside a cavity are usually used 
to generate squeezed light (e.g. [hi. [lol fl(| ) . 

In a waveguiding geometry that profits from a strong 
spatial confinement of the interacting optical fields in the 
transverse plane resulting in high values of the effective 
nonlinearity, another method to reach a nonlinear phase 
mismatch is possible. One of the nonlinearly interact- 
ing fields can be coupled through its evanescent waves 
into another field of the same frequency propagating in 
a neighbouring waveguide. An exchange of energy be- 
tween these two linearly coupled fields introduces a spa- 
tial modulation of the nonlinearly interacting field that 
can be set to compensate for the nonlinear phase mis- 
match [l7j . Interaction of fields in different waveguides 
through their evanescent waves can be used in various 
configurations that modify nonclassical properties of op- 
tical fields emerging from nonlinear interactions [l8|. 

A waveguiding geometry allows another possibility to 
tailor the nonlinear process - a linear periodic corruga- 
tion of the waveguide surface can be introduced as a 
distributed feedback that scatters the propagating fields 
19, 20]. Under suitable conditions, the scattered contri- 
butions interfere constructively and increase electric-field 
amplitudes of the propagating fields. This then results 
in higher conversion efficiencies of the nonlinear process 
[2lL I22I ] . Also squeezed- light generation is supported in 
this geometry, as discussed in [23| where scattering of the 
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second-harmonic field on the corrugation has been ne- 
glected. A waveguiding geometry with a periodic corru- 
gation can also be conveniently used for second-harmonic 
generation in Cerenkov configuration [24[. A general 
model of squeezed-light generation in nonlinear photonic 
structures has been developed in (25l |. 

In this paper, we show that scattering of the interact- 
ing fields caused by a linear periodic corrugation supports 
the generation of squeezed light in both the pump and 
second-subharmonic fields provided that an electric-field 
amplitude of at least one of the interacting fields is in- 
creased inside the waveguide as a consequence of scatter- 
ing. We note that an increase of electric-field amplitudes 
in the area of periodic corrugation is small compared with 
that occurring in layered photonic band-gap structures 
(with a deep grating) [H, [23, HM| - 

The generation of squeezed light has been in the cen- 
ter of attention in quantum optics for more than twenty 
years. It has been shown that the squeezed light can 
be generated also in Kerr media, nonlinear process of 
four-wave mixing, and diode lasers pumped by a sub- 
Poissonian current (for a review, see [HI). The level of 
squeezing observable in recent experiments [29L l30j ap- 
proaches 10 dB below the shot-noise level. 

The paper is organized as follows. In Sec. II, a quan- 
tum model of the nonlinear interaction including both 
Heisenberg equations for operator electric-field ampli- 
tudes and model of a generalized superposition of sig- 
nal and noise are presented. Conditions for an efficient 
squeezed-light generation are derived in Sec. III. A de- 
tailed analysis of the waveguide made of LiNbC>3 is con- 
tained in Sec. IV. Conclusions are drawn in Sec. V. Ap- 
pendix A is devoted to mode analysis of an anisotropic 
planar waveguide. 



x=z 



LiNt>03 




substrate z y LiNb03 

FIG. 1: Four optical fields interact in a nonlinear waveguide 
of thickness t and length L with a periodically-poled (pe- 
riod A„i) x' 2 ' susceptibility; A PF , A PB , A Sp , and A 3B mean 
forward-propagating pump, backward-propagating pump, 
forward-propagating second-subharmonic, and backward- 
propagating second-subharmonic electric-field amplitudes, re- 
spectively. A linear corrugation with depth ti and period A; 
is fabricated on the waveguide upper surface. Profile of the 
waveguide in the x — y plane is rectangular with width Ay 
and depth t; Ay 2> t is assumed. The waveguide is made of 
LiNbOs with the optical axis oriented along the x axis. 



the medium. The symbol A stands for Laplace opera- 
tor, (•) means a scalar product, and • denotes tensorial 
multiplication. Every spectral component of the element 
eij(r,o;) of the relative permittivity tensor in the consid- 
ered waveguide can be expressed as follows: 

Eij (x,y,z,uj) = lij (x, y, uj) [1 + Ae rj (x,y,z, uj)} . (2) 

Small variations of permittivity e described by 
Aeij(x, y, z, uj) are induced by a periodic corrugation of 
the waveguide surface. These variations of the elements 
Aeij(x, y, z, uj) along the z axis can be conveniently de- 
composed into harmonic waves: 



Aeij(x,y, z,ui) = e i j, q (x,y ) u))exp 



2tt ' 



(3) 



II. MODEL OF THE NONLINEAR 
INTERACTION 

An overall electric-field amplitude E(r, t) describing 
an optical field in the considered anisotropic nonlinear 
waveguide (shown in Fig. [1} is composed of two con- 
tributions; pump (or fundamental) electric-field ampli- 
tude E p (r, t) at frequency uj and second-subharmonic 
electric-field amplitude E s (r,i) at frequency w/2; i.e. 
E = E p + E s . We note that in case of second-harmonic 
generation, the field with the amplitude E p (r, t) is called 
second-harmonic and that with the amplitude E s (r,i) 
is known as pump. We keep the terminology used for 
second-subharmonic generation throughout the paper. 
The electric-field amplitude E obeys the wave equa- 
tion inside the waveguide with a nonlinear source term 
MM: 



AE- V(V • E) - /ie e 



d 2 E 



d 2 Pm 

1 



(1) 



In Eq. |T]) fi denotes vacuum permeability, eo vacuum 
permittivity, and P n i describes nonlinear polarization of 



where e q y are coefficients of the decomposition and A; 
is a period of the linear corrugation. Amplitude of the 
nonlinear polarization P n i of the medium is determined 
using tensor d of the second-order nonlinear coefficient: 



P nl (r,t)=2erjd(r)-E(r,i)E(r,i). 



(4) 



Taking into account geometry of our waveguide elements 
dyfe(r) of the nonlinear coefficient can be expressed as 
follows: 



d ijk (x,y,z) 



q— — oo 



d ijk ^(x,y)eyip 
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(5) 



where A„; describes the period of a possible periodical 
poling of the nonlinear material. 

The electric-field amplitudes of pump (E p ) and second- 
subharmonic (E s ) monochromatic waves can be ex- 
pressed in the form: 

E a (x,y,z,t) = i[A ap (z)e a (x,y)exp(i(3 a z-iuj a t) 

+ A aB (z)e a (x, y) exp(-i/? a z - iui a t) 



-H.c] 



(G) 
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where the symbols e p and e s refer to mode functions in 
the transverse plane of the beams. Amplitudes A pF and 
A SF [A Pb and A Sg ] describe forward- [backward-] propa- 
gating pump and second-subharmonic fields and are such 
that the quantities |v4 pj; .| 2 , |^4p B | 2 , |^s f | 2 j an d |A;b| 2 
give directly the number of photons in these fields. Sym- 
bol f3 a means propagation constant along the z axis in 
mode a whereas u a stands for frequency of this mode. 
Symbol H.c. replaces hermitian conjugated terms. We 
note that the considered anisotropic waveguide supports 
only TM guided modes; TE modes are not guided and 
thus do not contribute significantly to nonlinear interac- 
tion. 

Mode functions e p and e s describe the transverse pro- 
files of pump and second-subharmonic fields, respectively, 
and fulfill the following equations: 

V x (V x [e a (x,y)exp(±i/3 a z - iu a t)\) 
uo 2 

- -^eij(x,y, U! a ) ■ [e a (x, y) exp(±i(3 a z - ito a t)] = 0, 

a = p,s. (7) 

We note that mean values of permittivities (x, y, uj a ) 
are used for the determination of mode functions. The 
mode functions e p and e s are normalized to describe one 
photon in a mode inside the waveguide (see Appendix 
A). 

Substitution of Eqs. (JUH]) into Eq. ([T]) assuming 
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for a = p,s and b = F, B (analog 

of the slowly-varying envelope approximation to spatial 
evolution) results in the following equations for ampli- 
tudes A pF , A pB , A SF , and A SB : 

dA, w 



dz 

dA Sl 
dz 

dA P[ 
dz 

dA Pl 
dz 



iK s exp(—i8 s z)A SB 
+ 4K n i iq exp(i6 n i tq z)A PF A* F , 

-iK* exp(iS s z)A SF 

- 4K n i iq exp(-iS n i iq z)A PB A* E 

iK p e~xjp(—i5 p z)A PB 

- 2K* l q exp(-iS n i, q z)A 2 SF , 
—iK* exp(iS p z)A PF 

+ 2K nl,q exp(iS n l yq z)A 2 B . 



(8) 



Because the waveguide is made of an anisotropic ma- 
terial, also the following relation has to hold in order to 
derive correctly equations in Eq. (for details, see [3l1]): 



d 2 A ab (z) 
dz 2 

dAq b {z) 
dz 

, dA ab (z) 



dz 



e a (x,y)] z exp(±ip a z) z 
V ([e a (x, y)] z exp(±i/3 a z)) 
(V ■ [e a (x, y) exp(±if3 a z)]) z w 0, 



P,s, 



F.B, 



(9) 



and z stands for a unit vector along the z axis. We note 
that Eqs. © describe also a nonlinear coupler composed 
of two waveguides made of media whose modes in- 
teract through evanescent waves. Nonclassical properties 
of light propagating in this coupler have been studied in 
33L 3^ | . A scheme that allows to decompose interactions 
in this coupler into a sequence of fictitious interactions 
has been suggested in |3o| . 

Phase mismatches 5 S , 5 P , and 5 n i A occurring in Eq. ([8]) 
are given as follows: 



2tt 

S a = 2j3 a - — , a = p, s, 



5nl,q = f3 p -2(3 s +q 



2tt 



(10) 



Coefficient q equals ±1 for a periodically poled nonlinear 
material, whereas q = for a material without period- 
ical poling. Linear coupling constants K p and K s are 
determined along the expressions: 



u) 2 a J dxdyei(x,y,uj a ) ■ e* a (x,y)e a (x,y) 



2c 2 (3 a 
a = p, s. 



f dxdy\e a (x,y)\< 



(11) 



Similarly, the following expression can be found for non- 
linear coupling constants K n i A for q = 0, ±1: 



K„i, 



2iuj 2 J dxdyd q (x,y) ■ e p (x,y)e*(x,y)e*(x,y) 
C 2 I3 S f dxdy \e s (x,y)\ 2 

2iu 2 J dxdyd q (x,y) ■ e p (x,y)e* s (x,y)e* s (x,y) 



c 2 P P 



! dxdy\e p (x,y)\* 



(12) 



The last approximate equality in Eq. (|12p is valid pro- 
vided that to s /p s w ujpj '{3 P and due to the normalization 
of the mode functions e s and e p . This approximation as- 
sures, that only one nonlinear coupling constant occurs 
in Eqs. ([5]) which is important in quantum description. 

Expressions for linear (K s , K p ) and nonlinear (K n i y o, 
K n i,i) coupling constants appropriate for the considered 
waveguide and derived from Eqs. (fTTj) and (|12p can be 
found in Appendix A in Eqs. (|A10IIA12p . 

Quantum model of the nonlinear interaction in the con- 
sidered waveguide can be formulated changing the clas- 
sical envelope electric-field amplitudes A PF , A PB , A SF , 
and A Sg occurring in Eq. ([6]) into operators denoted as 
A PF , A PB , A SF , and A SB , respectively. A quantum ana- 
log of the classical equations written in Eq. ([S]) can then 
be derived from the Heisenberg equations (for details, see 



dX 

dz 



G, X 



(13) 



considering the following momentum operator G: 



G = 



TiK s exp(-i5 s z)Al A s 
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+hK p exp(-i5 p z)Al F A PB + H.c. 
- 2i?iK n i iq exp(i8 nl7q z)Al 2 F A pF 
+2ihK n i tq eKp(-i5 n i >q z)Al 2 B Ap B 



H.c. 



,(14) 



where q = or q — ±1. In Eq. (fl"3")l . fi means the reduced 
Planck constant, X stands for an arbitrary operator and 
the symbol [, ] denotes commutator. 

The operator quantum equations analogous to those 
written in Eq. ([5]) can be solved using the method of 
a small operator correction (denoted as 5 A) to a mean 
value (denoted as A) in which an operator electric-field 
amplitude A is decomposed as A = A + SA. This 
method provides a set of classical nonlinear equations for 
the mean values A that coincides with the set given in 
Eq. ([8]). The operator electric- field amplitude corrections 
8 A fulfill the following linear operator equations: 



dSA s 
dz 

dSA s 
dz 

dSA 



K S SA S 



K-F,q 



A PF SA SF 



A PB $A* S 



A* SF SA PF 



A* 5A„ 



PF 



dz 
dSA n 



K, p SA PB — K,p q A SF 8A Sp , 



dz 



1C*5A PL 



K* Bq A SB 5A s 



(15) 



The functions JC S , K p , and IC n i yq introduced in Eqs. JT 
are defined as: 



K. a = iK a exp(-iS a z), 
/C F , q = 4K n i iq ex.p{iS n i iq z), 
/Cb.q = 4:K n i iq ex.p(-i6 nl , q z). 



a = s,p, 



(16) 



Solution of the classical nonlinear equations written in 
Eqs. ([S]) can only be reached numerically using, e.g., a 
finite difference method called BVP [38]. This method 
requires an initial guess of the solution that can be con- 
veniently obtained when the nonlinear terms in Eqs. ([8]) 
are omitted. Then, the initial solution can be written as 
follows: 



4 0) 

ap 



A(0) 



exp —i 



.S a z 



exp 



.S a Z 



B n 



B n 



B a cos(A a z) + B a sin(A a 2) 



— ^- cos(A a z) + i— - sin(A a z; 
2K a K a 

~- sin(A a z) - i^- cos(A a z~ 
2K a K a 



s,P, 



and 



A„ 



4 



\K n 



s,p. 



(17) 



(18) 



In Eqs. (fT"7]) . constants B p , B p , B s , and B s are set ac- 
cording to the boundary conditions at both sides of the 
waveguide. 

We note that any solution of the nonlinear equations in 
Eqs. ((HI) obeys the following relation useful in a numerical 
computation: 



dz 



\A SF \ 2 + 2\A PF f-\A SB f-2\A PB \ 2 )=0. (19) 



The solution of the system of linear operator equations 
in Eqs. (| 1 5[) for the operator electric-field amplitude cor- 
rections 8 A can be found numerically and put into the 
following matrix form: 



SA.F,out 

SAbau 



( Uff Ufb \ 
\Ubf Ubb J 



where 



SAi 



Ms. in = 



fsA SF (o)\ 

sA PF (o) 
Uit F (o)/ 

fSA SB (L)\ 
SAUL) 



SA 



F.out 



5A PB {L) 
1< 



SA 



B,out 



\sAUl)J 



5AF,in 
SAB.out 



/SA SF (L)\ 
5At F (L) 
SA PF (L) 

y5Al(L)/ 

fSA SB (0)\ 
SAt B (0) 
SA PB (0) 

\<at B (o)/ 



(20) 



-(21) 



and L means the length of the waveguide. Matrices Uff, 
Ufb,Ubf, and Ubb are determined using the numerical 
solution of Eqs. (|15p . 

The following input-output relations among the oper- 
ator amplitude corrections SA, 



SAp^out 
SA B .out 



( Uff — UfbU b b Ubf UfbU b b 
V -U bb Ubf U bb 

SAf,i^ 
SAbm 

u ( 



(22) 
(23) 



are found solving Eqs. (|2T|) with respect to vectors 
SAF,out and S Ab, out- The output operator electric-field 
amplitude corrections contained in vectors SAF,out an d 
SAB,out obey boson commutation relations provided that 
the input operator electric-field amplitude corrections 
given in vectors <5AF,in and 5Ab,itx are ruled by boson 
commutation relations. We note that also certain com- 
mutation relations among the operator electric-field am- 
plitude corrections in vectors i5-4F,out and SAB,in can be 
derived (for details, see [37j). 

We restrict our considerations to states of optical fields 
that can be described using the generalized superposition 
of signal and noise 0. Thus coherent states, squeezed 
states as well as noise can be considered. Parameters 
Ba, Cj, D-jk, and ZXfc defined below are sufficient for the 
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description of any state of a two-mode optical field in this 
approximation [18j |: 



Overall phase-matching of the nonlinearly 
interacting fields 



D jk 



= (AA]AAj), 

= ((Ai,) 2 ), 

= (Ai,Ai fe ), 

= -(AitAi fc ) 



(24) 



AAj = Aj — (Aj); symbol ( ) denotes the quantum sta- 
tistical mean value. Expressions for the coefficients Bj, 
Cj, Djk, and Djk appropriate for outgoing helds can be 
derived (37j using elements of the matrix U given in Eq. 
([23]) and incident values of the coefficients Bj^A and 
Cj, m ,A related to anti-normal ordering of field operators 
(for details, see [H[): 



B j,in,A = cosh 2 (r j) +n chtj , 
Cj,in : A = ^cxp(z^ i )sinh(2r i ). 



(25) 



In Eq. (US]), Tj stands for a squeeze parameter of the 
incident j'-th field, "dj means a squeeze phase, and n c hj 
stands for a mean number of incident chaotic photons. 
Coefficients Djk t \ n ,A and Djk t \ n ,A f or an incident field are 
considered to be zero, i.e. the incident fields are assumed 
to be statistically independent. 

The maximum attainable value of squeezing of electric- 
field amplitude fluctuations is given by the value of a 
principal squeeze variance A 0. Both single- mode prin- 
cipal squeeze variances Xj and compound-mode principal 
squeeze variances \j (characterizing an overall field com- 
posed of two other fields) can be determined in terms of 
the coefficients Bj,Cj, Djk, and Djk given in Eq. ([24]) 
(for details, see, 



A, = l + 2[B J -\C j \], 



= 2[1 + 
-\Cj 



C, 



B k - 2^(D jk ) 
+ 2D jk \}; 



(26) 
(27) 



symbol means the real part of an expression. Values 
of the principal squeeze variance Aj (\jk) less than one 
(two) indicate squeezing in a single-mode (compound- 
mode) case. 



III. SUITABLE CONDITIONS FOR 
SQUEEZED-LIGHT GENERATION 

It occurs that there exist two conditions for an ef- 
ficient squeezed-light generation. The first condition 
comes from the requirement that the nonlinear interac- 
tion should be phase-matched along the whole waveguide, 
whereas the second one gives optimum conditions for the 
enhancement of electric-field amplitudes of the interact- 
ing optical fields inside the waveguide. 



Conditions for an optimum phase-matching of the in- 
teracting fields can be revealed, when we write the dif- 
ferential equation for the number N a of photons in field 
a; N a = A^A a . Using Eqs. ([8]) we arrive at the following 
differential equations: 



dN SI 
dz 



dNsi 
dz 



dN.. 



VF 



dz 

dN Pl 
dz 



-2%{K 3 eM-Msz)A* sp A SB } 

exp(iS n i, q z)A SF Ap F } , 
-2%{K s eM-M s z)A* p A SB } 

- 85i {K nhq exp(-iS n i, q z)A* s 2 B A PB } , 
-23 {K p exp(~iSpz)Ap F Ap g } 

- A^{K nhq e^{i5 nUq z)Af F A PF } , 
-23 {Kp e~K.p(~iS p z)A* F Ap B } 

+ m {K nhq exp{-i5 nliq z)Af B A PB } ; (28) 



symbol 3 denotes the imaginary part of an expression. 
We note that the first terms on the right-hand side of 
the first and the second (as well as the third and the 
fourth) equations in Eqs. (f2"5|) have the same sign because 
of counter-propagation of the fields [see also Eq. ([15])]. 
The nonlinear interaction described by the second terms 
on the right-hand sides of Eqs. (|28|) is weak and so we can 
judge the contribution of these terms using a perturba- 
tion approach. In the first step we neglect the nonlinear 
terms in Eqs. ([5]) and solve Eqs. ([5]) for field amplitudes 



A^J, Asb, A^J, and Ap"^. Then we insert this solution 



, (o) 



, (o) 



into the nonlinear terms in Eqs. (|28[) and find this way 
optimum conditions that maximize contributions of these 
terms. The solution for amplitudes A^Jt , As B , A P ° F , and 

Ap B coincides with that written in Eqs. (|17[) as an initial 
guess for the numerical solution and we rewrite it into 
the following suitable form: 



.5 a z 



A^(z) = exp -i^ 



x exp(iA a z) + B ap cxp(-i A Q z)] , 
B a — iB a 



B+ = 

2 

Bg + iBg 



Af B {z) = 



exp i 



.5 a z 



[B+ B exp(iA a z) + B aB exp(-iA a z)] 
5+ = - 6a .t 2Aa (B a -iB a ), 



4K n 



4K a 



-(B a + iB a ),a = p,s. 
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(29) 

Provided that a linear corrugation is missing in field a 
the solution for amplitudes and A„b can be obtained 
from the expressions in Eqs. (|29p using a sequence of two 
limits; S a — > 0, K a — > 0: 



B a 
B„ 



(30) 



The nonlinear interaction between the forward- 
propagating fields is described in Eqs. Ij28|) in our per- 
turbation approach by the term 

[K nhq exp(z<5 n ^z)4 / 2 4° F ) } (31) 

that, after substituting the expressions for amplitudes 
A f) and A\ 
eight terms: 



A^sp and Ap} from Eqs. ([29]) . splits into the following 



x exp(i[£„ i>9 + <5 S - 5 p /2 t A s tA s ± A p ]z) L(32) 



The nonlinear interaction is efficient under the condition 
that one of these terms does not oscillate along the z 
axis. This gives us eight possible conditions that com- 
bine nonlinear phase-mismatch 8 n i^ q and parameters of 
the corrugation 5 S , S p , K s , and K p : 



5nl,q 



S p /2 t A s t A s ± A p 



0. 



(33) 



It depends on a given waveguide and initial conditions 
which out of these eight conditions leads to an efficient 
nonlinear interaction. We note that the conditions in 
Eq. (|33p are valid also for the nonlinear interaction be- 
tween the backward-propagating fields characterized by 

the term $1 (K„i, g exp(-i8 nhq z)A ( S A p }} in Eqs. (J2SJ). 

We consider two special cases in which a periodic cor- 
rugation is present either in the pump or the second- 
subharmonic field. Assuming the corrugation in the 
pump field the conditions in Eq. (|33[) get the form: 



$ni, q - Sp/2 ± A p = 0. 



(34) 



Sign - (+) is suitable for 8 n ^ q > {5 n i, q < 0) when we 
solve Eq. ((311) for 8 p : 



Sp — $nl,q 



\K P \ 2 . 



8 n . 



l,q 



(35) 



he. 8 n l,q an< i S p have the same sign. The expression in 
Eq. (|3"B")) then determines the period A; of linear corruga- 
tion: 



A, 



/3 p -S n i >q /2-\Kp\y(28 nltq )- 



(36) 



On the other hand, the conditions 

8 nUq + S s /2 T A s T A s = 



(37) 



are suitable for the periodic corrugation in the second- 
subharmonic field. When S n i >q > (S n ^ q < 0) signs - (+) 
in Eq. (|37p are appropriate and we have: 



6 S 



*nl,q 



2\K„ 



0nl,q 



(38) 



i.e. 5 n i tq and S s have the opposite sign. The period A/ of 
linear corrugation is then given as 



A, 



7T 



Ps + Snl,q/4+ \K P \ 2 /S n i iq ' 



(39) 



B. Enhancement of amplitudes of the interacting 
fields 



The greatest increase of electric-field amplitudes in- 
side the waveguide occurs under the condition of trans- 
parency for the incoming field [2lj]. We note that this 
property also occurs in layered structures with band-gaps 
and transmission peaks where the field is well localized 
inside the structure provided that it lies in a transmission 
peak [II]. 

A transmission peak of the waveguide can be found 
from the condition that a backward-propagating field is 
zero both at the beginning and at the end of the waveg- 
uide, because it is not initially seeded. These require- 
ments are fulfilled by the solution in Eqs. (fl~7|) under the 
following conditions: 

mn 

A a = — , to = 1,2,...; 

B a = A^(0); B a = ^-B a , a=p 7 s. (40) 

A natural number m counts areas of transmission. 

The linear phase mismatch S a is determined from the 
first equation in Eqs. (f40|) in the following form 



5 a = ±2 



yJ(^y + \K a \*, a = p,s, (41) 



and the corresponding period A; of linear corrugation is 
given as: 



A, 



(3 a ± y/(mir/L)2 + \K a \< 



(42) 



The conditions in Eqs. (|35|) and (|41j) for a periodic 
corrugation in the pump field can be combined together 
to provide the following formula for the coupling constant 



(43) 
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The linear phase mismatch S p is then determined along 
the formula in Eq. (|4"Tj) such that the signs of S n ^ g and 
dp are the same. 

Similarly the conditions for a periodic corrugation in 
the second-subharmonic field written in Eqs. (|38p and 
(141) and considered together lead to an expression for 
the coupling constant K s : 

The sign of linear phase mismatch 5 S determined from 
Eq. (|4Pj) is opposite to that of 5 n i, q - 

IV. SQUEEZED-LIGHT GENERATION - 
NUMERICAL ANALYSIS 

The discussion of squeezed-light generation is decom- 
posed into three parts. In the first one, we pay a 
detailed attention to attainable characteristics of the 
considered waveguide. The second part is devoted to 
second-subharmonic generation, i.e. an incident strong 
pump field is assumed. In the third part, a strong in- 
cident second-subharmonic field is assumed, i.e. second- 
harmonic generation is studied. 



A. Characteristic parameters of the waveguide 

We consider a waveguide made of LiNbC>3 with length 
L = 1 x m and width Ay = 1 x 1CP 5 m pumped 

at the wavelength \ p = 0.534 x 1(T 6 m (A s = 1.064 x 
10~ 6 m). More details are contained in Appendix A. 
We require a single-mode operation at both the pump- 
and second-subharmonic-field frequencies that can be 
achieved for values of the thickness t of the waveguide 
in the range t G (0.403, 0.544) x 10~ 6 m. We also ne- 
glect losses in the waveguide that may be caused both by 
absorption inside the waveguide and scattering of light 
on the corrugation that does not propagate into guided 
modes. In practise, material absorption is weak at the 
used wavelengths. Scattering of light on imperfections 
of both the guided structure and periodic corrugation 
determines losses in a real waveguide and leads to degra- 
dation of squeezing. We also assume that intensity of the 
pump field is such that self-phase modulation due to x 3 
nonlinearity does not occur. 

Values of the natural nonlinear phase mismatch 8 n i t o 
are high (around 1.7 x 10 6 m _1 ) in the region of single- 
mode operation and depend on the thickness t of the 
waveguide (see Fig. [2^,). An additional periodical poling 
has to be introduced to compensate, at least partially, 
for this mismatch and allow an efficient nonlinear pro- 
cess [23j]. Values of the nonlinear coupling constant Kj^ja 
depicted in Fig. increase with the increasing values of 
the thickness t because of the thicker the waveguide the 
more the modes are localized inside the waveguide and 
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FIG. 2: Natural nonlinear phase mismatch 8 n i,o (a) and ab- 
solute value of the nonlinear coupling constant K n i$ (b) as 
they depend on the thickness t of the waveguide are shown in 
the region of single-mode operation; constant K n l,o is deter- 
mined for amplitudes that correspond to one pump and one 
second-subharmonic photon inside the waveguide. 

so the greater the overlap of the nonlinearly interacting 
electric-field amplitudes. Attainable values of the lin- 
ear coupling constants K p and K s for different values of 
thickness t and depth ti of corrugation can be obtained 
from Fig. [3l Values of the coupling constants K p and 
K s for a given thickness t increase monotonously with 
the increasing values of the depth ti of corrugation. For 
small values of the thickness i, values of the coupling 
constant K s are small because the waveguide is thin for 
the second-subharmonic field and so a considerable part 
of the field is outside the waveguide and cannot be scat- 
tered by the corrugation. 

Amplitudes \A ah \ of the dimensionless incident strong 
electric-field amplitudes are determined from the incident 
power P ab along the relation: 

l^»l = y^yS a =P> s > b = F,B. (45) 

On the other hand, the power P° b ut of an outgoing field 
is given as follows: 

hi . 2 

pout _ nuJ a I a |2 

ab ' p a L 1 abl 

huj 2 

= lTfN ab , a = p,s, b = F,B; (46) 

N ab denotes the number of photons leaving the waveg- 
uide. 

New dimensionless parameters are convenient for the 
discussion of behavior of the waveguide; 

z r = z/L, 
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|K | (10W 1 ) 
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b) 
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FIG. 3: Contour plot of the absolute value of pump-field 
[second-subharmonic-field] linear coupling constant K v [K s ] 
as it depends on thickness t of the waveguide and depth ti 
of linear corrugation is shown in (a) [(b)] in the region of 
single-mode operation. 



Af = Ai/L,A r nU = A n ijL, 
P: = L(3 a ,6 r a = LS a ,K: = LK a , 



°nl,q 



LS. 



LK 



nl,q- 



(47) 



Applying these parameters the waveguide extends from 
z r = to z r — 1. The dimensionless parameters enable 
to understand the behavior of the waveguide as it de- 
pends on the length L using the graphs and discussion 
bellow. 



B. Second-subharmonic generation 

As a reference for the efficiency of squeezed-light gen- 
eration we consider the waveguide with periodical poling 
and assume that it is pumped by the power of 2 W. The 




3.52 3.54 3 56 3.58 

10 3 A„, r 



FIG. 4: Principal squeeze variances X SF (solid curve) and 
X PF (solid curve with *) as functions of the period A 1 nl of 
poling for second-subharmonic generation; to understand the 
behavior with respect to the length L, the period A n i should 
be decomposed as A°, + AA„;, then S 7 nl q « — 2TvqAA 7 nl /A^ ll 2 , 
A° nl = 2nq/{2p s - f3 p ), AA r nl = LAA nl ; i = 5x 1(T 7 m, 
K p = 0, K T S = 0, P PF = 2 W, P SF = 1 x 1(T 10 W (a negligible 
seeding to substitute the spontaneous process in the classical 
equations), arg(A PF ) = 0, a.rg(A SF ) = 0, A p 



0, A B 



0; 



A" 



3.552 x 10" 6 m 



nonlinear interaction is perfectly phase matched for the 
period of poling A r nl w 3.552 x 10~ 3 where we have for 
the principal squeeze variances A SF ~ 0.45 and A PF ~ 1 
(see Fig. |4]). The more distant the value of A r nl from the 
above-mentioned optimum value is, the larger the non- 
linear phase mismatch and the larger the value of the 
principal squeeze variance \ SF . 

Considering a periodic corrugation in the pump field 
with such parameters that enhancement of the pump field 
inside the waveguide occurs, better values of squeezing in 
the second-subharmonic field can be reached. However, 
nonzero values of the nonlinear phase mismatch 8 n i t \ are 
important to reach better squeezing because they have 
to compensate periodic spatial oscillations caused by the 
corrugation (see [1(3]). A perfect phase matching of all 
the processes occurring in the waveguide can be reached 
this way [see Eq. ([53])]. 

A typical dependence of the principal squeeze variance 
X SF as well as the number N SF of photons leaving the 
waveguide for forward-propagating second-subharmonic 
field for attainable values of parameters of the corruga- 
tion is shown in Fig. [5] assuming a fixed value of the 
nonlinear phase mismatch 8 r nl 1 equal to - 10.82 {A r nl = 
3.53 x 10~ 3 ), i.e. it corresponds roughly to the second 
local maximum of A SF in the curve in Fig. 2] We can 
clearly see that an efficient nonlinear interaction occurs in 
strips that correspond to transmission peaks; the larger 
the number m of a transmission peak [see Eq. (|40[)] the 
weaker the effective nonlinear interaction. The neces- 
sity to fulfill also the condition for perfect phase match- 
ing given in Eq. (|35[) is evident. The principal squeeze 



variance \ SF reaches values around 0.2 inside the strips 
around the first several transmission peaks. The larger 
the number to of a transmission peak, the greater values 
of linear coupling constant and linear phase mismatch 
5p have to be used to reach high levels of squeezing. Up to 
several forward-propagating photons can be present in- 
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FIG. 5: Contour plots of principal squeeze variance X SF (a) 
and number N SF of photons leaving the waveguide (b) for 
forward-propagating second-subharmonic field as they depend 
on parameters \Kp\ and 5p of the corrugation in the pump 
field; more-less equidistant curves in the graphs indicate po- 
sitions of transmission peaks [Eq. (|40|l ]. the last curve going 
up is given by the condition in Eq. (|35[1 ; h-ni = 3.53 x I CP 3 
(finl l = —10.82), arg(.ftTp) = it/2; values of the other pa- 
rameters are the same as in Fig. [4] Triangles in lower right 
corners with formally zero values in both graphs lie in the 
region with an exponential behavior of classical amplitudes 
that is not suitable for nonclassical-light generation. 

side the waveguide (see Fig. [SJd) at a given time instant. 
This means that the power of the outgoing field is of the 
order of 1CP 8 W (energy of one second-subharmonic pho- 
ton inside the waveguide of thickness t — 5 x 10~ 7 m and 
length L = 1 x 10~ 3 m corresponds to the output power 
of 2.58 x 10~ 8 W). Only the first and the second trans- 
mission peaks can give reasonable values of the power of 
the outgoing field. 

Effect of a periodic corrugation to the enhancement 
of the nonlinear interaction can be approximately quan- 
tified as follows. Because the pump field is strong, its 
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FIG. 6: Enhancement factor M as a function of linear 
coupling constant K r p for the first five transmission peaks: 
m — 1 (solid curve), 2 (solid curve with o), 3 (solid curve 
with *), 4 (solid curve with A), and 5 (solid curve with o); 
t = 5 x 1(T 7 m. 



depletion by the nonlinear process can be omitted and 
its amplitude A PF (z) along the waveguide is given by the 
formula in Eq. (|29[) . In our case, sign + in the condi- 
tion in Eq. (|34[) is valid and this means that the nonlin- 
ear process exploits efficiently the first term in Eq. ([29]) 
that is multiplied by the coefficient B+ . An effective en- 
hancement of the pump-field amplitude A PF inside the 
waveguide can be given by the ratio B pF /A PF (0) that we 
call an enhancement factor M.. The enhancement factor 
M. can be expressed as follows provided that the pump 
field lies in a transmission peak: 

2A p + <5 p 11 l \K p \2 + (mn/Ly 
M -—[K p —-2 + 2y WL)2 ' (48) 

According to Eq. ([48]) the greatest value of enhance- 
ment factor M. occurs at the first transmission peak 
(to = 1) and the greater the linear coupling constant K p 
the greater the value of enhancement factor M. . The en- 
hancement factor M. for the first five transmission peaks 
is shown in Fig. [6] as a function of the linear coupling 
constant K p . The greatest enhancement of an electric- 
field amplitude in the first transmission peak is accom- 
panied by the greatest difference between the maximum 
and minimum values of electric-field intensities along the 
waveguide. This poses the question whether other types 
of distributed feedback resonators (like a quarter wave 
shifted distributed feedback resonator) giving a more uni- 
form distribution of electric-field intensities along the 
waveguide j4l| can lead even to a better enhancement 
of the nonlinear process. Modelling of such distributed 
feedback resonators is, however, beyond the scope of the 
developed quantum consistent model and so we keep this 
question open. 

It is useful to compare the achievable values of the 
enhancement factor M with those typical for cavity ge- 
ometries. Assuming a symmetric planar cavity with 
mirrors having intensity reflection R, an electric-field 
amplitude A cav inside a cavity is determined from an 
incident electric-field amplitude A ml along the formula 
A cav = A lm / y/l — R. High quality cavities can have 
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R = 99.9 % and so the amplitude A ca,v is enhanced 
with respect to the amplitude A ml by factor 30, i.e. the 
enhancement of electric-field amplitudes is considerably 
greater in this case. On the other hand, the effective 
nonlincarity increases also due to field confinement in the 
transverse plane in a waveguide. 

Optimum values of waveguide parameters with respect 
to squeezed-light generation can be determined along the 
following procedure: 

• The greatest possible value of depth t\ of a peri- 
odic corrugation should be used to maximize the 
pump-field scattering. This gives the value of lin- 
ear coupling constant K p . From practical point of 
view, the depth t; of periodic corrugation is limited 
by technological reasons (also validity of the model 
would have to be judged for deeper corrugations). 

• The value of linear phase mismatch SL should be 
determined along the relation in Eq. (j41j) to insure 
that the pump field is in the first transmission peak. 
Both signs of the linear phase mismatch S p are pos- 
sible and the corresponding period A; of linear cor- 
rugation is determined according to Eq. (fT0| . 

• Solution of Eq. (|43p gives an appropriate value of 
the nonlinear phase mismatch S n i tq under the fol- 
lowing requirements: sign + in Eq. (|43[) is used and 
sign of the determined nonlinear phase mismatch 
5 n i,q has to be the same as that of the linear phase 
mismatch 8 P determined in the previous step. 

• Numerical analysis in the surroundings of the 
analytically-found values of parameters S p and 
S n i,±i finally gives the values of waveguide param- 
eters optimum for squeezed-light generation. 

This procedure is documented in Fig. [7J The depen- 
dence of the optimum value of nonlinear phase mismatch 
S^i i on the linear coupling constant K r v is shown in 
Fig. [Jji. The greater the value of linear coupling constant 
K p the greater the optimum value of nonlinear phase 
mismatch 5^ . The corresponding values of the princi- 
pal squeeze variance X SF are depicted in Fig. [JJs. The 
larger the value of nonlinear phase mismatch 8^ 1 the 
better the values of the principal squeeze variance A SF . 
Curves in Fig. [TJj also indicate that sign + in Eq. (j4"3"|) 
is appropriate and the first transmission peak gives the 
best values of principal squeeze variance X PF . Keeping 
the incident pump-field power fixed, the depth t\ of peri- 
odic corrugation limits the achievable values of principal 
squeeze variance X SF . 

To see usefulness of the corrugation we compare two 
configurations: a perfectly quasi-phase matched waveg- 
uide and a non-perfectly quasi-phase matched waveg- 
uide with a suitable corrugation that compensates for 
the given phase mismatch. The waveguide with corruga- 
tion gives better values of the principal squeeze variance 
A SF and also considerably greater values of the number 
N SF of photons leaving the waveguide, as documented 
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FIG. 7: Nonlinear phase mismatch 5£ ( x as a function of 
the absolute value \KZ\ of linear coupling constant [Eq. (143 |l ] 
(a) and principal squeeze variance X 3F as it depends on the 
nonlinear phase mismatch S , nl l (b); first (m=l) and second 
(m=2) transmission peaks as well as signs + and - in Eq. (|43[) 
are considered [m = 1+ (solid curve with *), m = 1— (solid 
curve with o), m = 2+ (solid curve with o), and m = 2— 
(solid curve with A)]; axg(Kp) = n/2; values of the other 
parameters are the same as in Fig. [4] 

in Fig. [5J Also improvement caused by an introduction 
of corrugation into a non-perfectly quasi-phase-matched 
waveguide is worth mentioning (see Fig. [5]). 

Benefit of the corrugation to squeezed-light generation 
can be quantified defining coefficient T> that gives the ra- 
tio (in dB) of the principal squeeze variance X SF reached 
with a corrugation and the principal squeeze variance X T SF 
characterizing a perfectly quasi-phase-matched waveg- 
uide without corrugation: 

P = -lOlog lo 0pf); (49) 

log 10 stands for decimal logarithm. The value of coeffi- 
cient T> determined under the optimum values of param- 
eters of the corrugation as given in Sec. Ill increases with 
an increasing incident pump- field power P pF (see Fig.[§]). 
According to curves in Fig. [5] the use of periodic corru- 
gation leads to a significant improvement of squeezing, 
especially for deeper corrugations and greater incident 
pump powers. 

The role of length L of the waveguide can be addressed 
using the above-mentioned results as follows. Because 
the nonlinear phase mismatch 8 n i,i does not depend on 
the length L, the dimcnsionless nonlinear phase mis- 
match Sj U 1 depicted in Fig. [7^, is linearly proportional 
to the length L and, according to curves in Fig. [7p, the 
larger the value of length L the lower the value of the 
principal squeeze variance X SF . The analysis of condi- 
tions in Eqs. (|4"Tj) and (|4"3")l for optimum squeezed-light 
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FIG. 8: Principal squeeze variance X SF (a) and number N SF 
of photons leaving the waveguide for forward-propagating 
second-subharmonic field as they depend on the incident 
pump-field power P PF under different conditions: perfect 
quasi-phase matching (<5£; 1 — 0) without corrugation (solid 
curve), quasi-phase matching with = 3.53 x 10 -3 (5„i i = 
— 10.82) without corrugation (solid curve with o), and quasi- 
phase matching with K r nl = 3.53 x 10~ 3 together with a corru- 
gation [its parameters are given by the conditions in Eqs. (|43p 
(m = 1) and (|35|) ] (solid curve with o); arg(A"p) = 7r/2; values 
of the other parameters are the same as in Fig. U 
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FIG. 9: Coefficient T> giving the improvement of squeezing 
as it depends on the incident pump power P PF for \K p \ — 5 
(solid curve with o), \K p \ = 10 (solid curve with *), \K P \ — 20 
(solid curve with A), and \Kp\ = 50 (solid curve with o). Solid 
curve without symbols is drawn for i = —10.82. Tuning to 
the first transmission peak (m = 1) as well as the condition 
for overall phase matching are assumed; t = 5 x 10 -7 m, 
K r s = 0, arg(^) = tt/2, P sf = 1 x 10" 10 W, arg(A aj? ) = 0, 
argfApp) =0, A SB — A PB =0. 



generation shows that S n i t i — > \K P \ and S p — > 2\K P \ in 
the limit of large length L. We note that the distance 
between the curves corresponding to adjacent transmis- 
sion peaks in contour plots depicted in Fig. [5] behaves as 
l/L. 

A periodic corrugation can be alternatively introduced 
into the second-subharmonic field. The results obtained 



for principal squeeze variances A and numbers N of pho- 
tons are comparable to those achieved with a corruga- 
tion in the pump field owing to an increase of electric- 
field amplitudes of the second-subharmonic field inside 
the waveguide. Suitable conditions are given in Eqs. (|38|) 
and (|40p in this case. 



C. Second- harmonic generation 

Second-harmonic generation requires a strong incident 
second-subharmonic field. As a reference we consider an 
incident power P SF of the second-subharmonic field to 
be 2 W and periodical poling giving perfect quasi-phase 
matching. We then have for the principal squeeze vari- 
ances A SF « 0.87 and X PF » 0.99 (A^ « 3.552 x 10~ 3 ). 

The nonlinearly interacting fields behave similarly as in 
the case of second-subharmonic generation. If we intro- 
duce a periodic corrugation into the second-subharmonic 
field and set the nonlinear phase mismatch 6^ l equal 
to 3.53 x 10 -3 , values of the principal squeeze variance 
A SF approach 0.6 inside the strips in the plane spanned by 
variables \K r s \ and 8 r s where the conditions from Eqs. [)38p 
and (f^D|) are fulfilled (compare Fig. [5]). The number 
N SF of second-subharmonic photons at the output of the 
waveguide (together with the output power) decreases 
by an order of magnitude inside these strips compared 
to the incident power as a consequence of transfer of en- 
ergy from the strong second-subharmonic field into the 
pump field (due to an efficient nonlinear interaction) and 
also transfer of energy into the backward-propagating 
field (due to scattering) is considerable. Despite this the 
squeezed second-subharmonic field remains very strong, 
it contains about 10 6 photons inside the waveguide. The 
pump field that is only weakly squeezed for perfect quasi- 
phase matching can reach values of the principal squeeze 
variance X PF around 0.8 assuming a corrugation with pa- 
rameters given by Eqs. (|3"8")l and (|40[) . The pump field 
gets a considerable amount of energy from the second- 
subharmonic field and so typical values of the number 
N PF of pump photons leaving the waveguide can reach 
10 6 ; i.e. the output power is of the order of 10 _1 W (en- 
ergy of one pump photon inside the waveguide of thick- 
ness t = 5 x 10~ 7 m and length L = 1 X 10~ 3 m corre- 
sponds to the output power of 5.13 x 10~ 8 W). 

A periodic corrugation can also be put into the pump 
field and we have qualitatively the same results as if the 
corrugation is present in the second-subharmonic field. 
Assuming i = —10.82 as above, values of the principal 
squeeze variance X SF can reach even 0.3. On the other 
hand, values of the principal squeeze variance X PF lie 
above 0.9. 

The influence of phases of the interacting optical fields 
to the nonlinear process is of interest. It can be easily 
shown that any solution of Eqs. (JSJ) depends only on the 
phase if) = SLTg(K p ) — 2 &rg(K s ). In our numerical inves- 
tigations, we did not observe any dependence of numbers 
N of photons as well as principle squeeze variances A on 
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the phase ip. On the other hand, these quantities depend 
weakly on phases arg(A SF ) and &rg(A pF ) of the incident 
fields. However, this dependence is very weak under the 
conditions where a strongly squeezed light is generated. 

At the end, we compare values of squeezing achievable 
in the considered waveguide with those reached in the 
commonly used cavity geometries. The best achieved val- 
ues of squeezing generated in a cavity approach to 10 dB 
[30| below the shot-noise level at present, i.e. values of 
the principal squeeze variance A SF lie slightly above 0.1. 
The analyzed waveguide does not reach so low values of 
the principal squeeze variance A SF (X SF « 0.2) because, 
as pointed out above, the enhancement of the pump-field 
amplitude inside the nonlinear medium is considerably 
lower compared to high quality cavities. On the other 
hand, the considered waveguide is relatively broad along 
the y axis (Ay = 1 x 10 -5 m) and so narrowing of the 
waveguide is possible. This would lead to greater values 
of the effective nonlinearity and subsequently to better 
values of squeezing. Also longer waveguides can be con- 
sidered. We believe that the analyzed waveguide with 
a periodic corrugation has potential to deliver squeezed 
light with values of parameters comparable to those mea- 
sured in cavity geometries. 



V. CONCLUSIONS 

We have shown that an additional scattering of two 
nonlinearly interacting optical fields caused by a small 
periodic corrugation on the surface of the waveguide can 
lead to an enhancement of the nonlinear process thus 
resulting in higher generation rates and better values of 
squeezing. Origin of this enhancement lies in constructive 
interference of the scattered fields leading to higher values 
of electric-field amplitudes inside the waveguide. Opti- 
mum conditions for this enhancement have been found 
approximately analytically and confirmed numerically. 
To observe this effect the natural nonlinear phase mis- 
match has to be nonzero in order to match with periodic 
oscillations caused by scattering at the corrugation. The 
deeper the corrugation and the higher the incident pump 
power the lower the values of the principal squeeze vari- 
ances. A periodic corrugation can be designed to match 
either the pump or the second-subharmonic field, or even 
both of them. 

The obtained results have shown that nonlinear planar 
waveguides with a periodically corrugated surface repre- 
sent a promising source of squeezed light for integrated 
optoelectronics of near future. 



APPENDIX A: MODES OF AN ANISOTROPIC 
WAVEGUIDE 

A mode of the considered waveguide [3£| depicted in 
Fig. [U is given as a solution of the wave equation written 
in Eq. ((7|). The waveguide is made of LiNb03 crystal us- 



ing the method of proton exchange. The crystallographic 
z axis coincides with the x axis of the coordinate system 
(see Fig.[I|. Ordinary (n Sj0 ) and extraordinary (n SiC ) in- 
dices of refraction of LiNbOa valid for the substrate and 
used in calculations are given as: 



B a 



D a X 2 , 



a 2 -a 

A = 4.91300, B = 0.118717, 
Co = 0.045932, D = 0.0278, 
A e = 4.57906,^ = 0.099318, 
C e = 0.042286, D e = 0.0224; 



a = o.e, 



(Al) 



wavelength A is in fim. After proton exchange, ordinary 
(n w _ ) and extraordinary (n w ^ e ) indices of refraction of 
LiNbOs characterizing the waveguide are reached: 



Tlw,o 
x2 



n 8 

Us 



(SnY = Ai 



1, 
Bi 



£iA 2 ; 



A 2 -Ci 

A 1 = 0.007596, £i = 0.001129, 

Ci = 0.116926, Di = -0.0003126. (A2) 

We assume that air is present above the waveguide, i.e.: 

n u = 1. (A3) 

Because only the extraordinary index of refraction 
of LiNbOs increases during proton exchange, only TM 
waves can be guided. For this reason, instead of solving 
Eq. ([7]) for x and z components of the electric-field mode 
functions e a , we solve the following equation for the only 
nonzero y component of the magnetic-field mode func- 
tions h (x) (fields are assumed to be homogeneous along 
the y axis) [311 ]: 



d 2 h a (x) 
dx 2 



e zz (x,uj a ) 2 

— / \ Pa 

e xx (x,uj a ) 



e zz (x,uj a )ujl 



K(x) =0. 

(A4) 



The solution of Eq. (|A4j) for y component of the 
magnetic-field mode function h a (ir) (a = p,s) can be 
written as: 



-C a ^r- exp(— qx), x > 0; 

la 



= c a 



•— cos(h a x) + sin(h a x) 
h a 



= -C a 

x exp(p a x), x < —t; 



— cos(h a t) + sin(/i Q i) 



, < x < -t; 
exp(p a t) 

(A5) 



where C a denotes a normalization constant. We have 
the following expressions for the coefficients h a , g a , p a , 
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p a and q a for the considered orientation of LiNbOs 



la 



Pa 



Pa 



<la 



(uj a )u> 

a \ f Tlw.o 



n s , (u>a) 

1 \> 



Pa) - 



n s , (uJ a )uJ a 



n w,o(Ug) 
n l,,o{Ua) 



Pa, 



-Qa- 



(A6) 



The solution written in Eq. (|A5j) holds provided that the 
following dispersion relation giving a propagation con- 
stant P a as a function of frequency uj a is fulfilled: 



t&n(h a t) 



hajpg + q a ) 
K ~ PaQa 



(A7) 



We note that possible solutions of Eq. (|A7|) for (3 a he in 
the interval n Ste (u} a )u)a/c < f3 a < n Wte {w a )w a /c. 

Components of the electric-field mode function e a (x) 
can be derived from the magnetic-field mode functions 
h a (x) along the relations: 



[e a (x)] x 
[e a (x)] y 
[e a (x)]z = 



Pa 



^a^-O^xx 



0. 



[K(x)]y, 



d[h a {3 



LUa£0£z Z (x,UJa) 



dx 



(A8) 



The normalization constants C a occurring in Eqs. (|A5[) 
are determined from the condition that the mode func- 
tions e a (x) describe one photon with energy huj a inside 
the waveguide (of length L and thickness Ay) : 



we have e±i — i(e— l)/(ire) in this case using Eqs. ([2]) 
and d3|). If the waveguide is periodically poled d±i = 
—2i/nd in Eq. ([5]) and the remaining coefficients may 
be omitted. Using the electric-field mode functions e p 
and e s determined in Eqs. (|A8[) . linear (K s , K p ) and 
nonlinear (K n i t o, K n i ; {) coupling coefficients defined in 
Eqs. (fTT|) and (fT2|) can be rearranged into the form: 



K a = 



2irc 2 f3 a 



n w .e.(.Ua) - 1 f° , |r / si ,2 

dx \ \& a {x)\ x \ 



n l,e(^a) J-t 



+ „2 (,,\ / dx \[ e a(x)\z\ 
rl w,o\ UJ a) J-t, 



dx(\[e a (x)} x \ 2 + \[e a (x)} z \ 2 ) 



a = p,s, 

■ 1 r0 

K„i,o = ^j- J dxd-e p {x)e*(x)e*(x) 

OO 



(A10) 



dx(\[e s (x)] x \ 2 + \[e s (x)] z \ 2 ) 



K„ii — 



21 K 

7T 



(All) 
(A12) 



Nonzero coefficients of the nonlinear tensor d of 
LiNbOs used in calculations are the following: 

3.1 x Kr 12 mV~\ 



41.05 x 10- 12 mV _1 



yyx = 5.87 x l(T 12 mV~\ 

(A13) 
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2e AyL dx e xx (x,u a ) |[e (x)] x |^ 



+ e zz (x,uj a ) \[e a (x)] z \" 



(A9) 



The corrugation on the surface causes periodic changes 
of values of permittivity e for x e (0, — tj) (see Fig.rjJ and 
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